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Abstract 

The contribution of the two-step process ~/d — > pn — > TT°d to the imaginary 
part of the amplitude for coherent pion production on the deuteron is calcu- 
lated exploiting unitarity constraints. The result shows that this absorptive 
process is not negligible and has to be considered in an extraction of the el- 
ementary neutron production amplitude from the jd — ► ir°d cross section at 
threshold. In addition, it is argued that a consistent calculation of 'yd — ► Tr°d 
in baryon chiral perturbation theory beyond next-to-leading order requires 

the inclusion of this absorptive process. 
PACS numbers: 21.45. +v, 25.20.Lj 
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Recently there has been considerable interest in the coherent electromagnetic production 
of pions from the deuteron near threshold. The main motivation thereby is to gain informa- 
tion on the elementary neutron amplitude jn — > 7T°n which is experimentally not directly 
accessible. A first measurement of the ed — > e7r°d reaction near q 2 = -0.075 GeV 2 /c 2 will be 
performed soon at MAMI |TJ. However, it is already known for a long time that at threshold 
the 7t£ — > ir°d process is dominated by two-nucleon production mechanisms [[|. Therefore, 
a careful theoretical analysis which allows the separation of the one-nucleon process is es- 
sential. Recently, chiral perturbation theory (xPT) in the heavy baryon formulation has 
been applied to this reaction by Beane et al. predicting a real threshold amplitude. Very 
recently, this work has been improved and extended beyond next-to-leading order in the 
chiral power counting scheme |4j]. 

It is the purpose of the present paper to point out that an additional contribution of an 
absorptive two-step process, 'yd — > pn — *> n°d, has to be included, which leads to a complex 
amplitude even at threshold. The presence of such a competing deuteron disintegration 
channel has an analogue in the case of nd elastic scattering, where the contribution of the 
absorptive process ird — > NN — > red is known to be of the order of 10% of the total amplitude 
(see e.g. 0). An effect of this size would not be negligible for the electromagnetic reaction 
due to the relative smallness of the single-nucleon amplitude, one is mainly interested in. In 
the case of nd — > ird, the imaginary part of the scattering length a n d is related to the total 
absorption cross section through 

Qma 7Td =^- lim p %d a(ird -> X), (1) 

where p n d is the pion momentum in the cm. system. This relation follows directly from 
the optical theorem. An analogous unitarity constraint on the 'yd — > ir d amplitude near 
threshold is present and will be treated as first point below. In view of the fact that the 
absorptive process has not been considered in [0,|J, we will analyze its role within the %PT 
framework as second point. 

As is well known, the unitarity of the S'-matrix leads to constraints for the corresponding 
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reaction amplitudes. For our purposes, it is sufficient to consider coupled two-particle chan- 
nels where the particles are subject to interactions which are invariant under time-reversal. 
In this case, following the conventions of Ref. ||, the imaginary part of the partial wave 
T-matrix element of total angular momentum J for the reaction a = (a^) —>■ b = (6162) 
fulfils the following unitarity constraint 

Here, A Ci is the helicity of particle i in the channel c, and p c denotes the cm. momentum 
of the two particles in the channel c. The first sum on the rhs of (|2|) runs over all open 
channels c for a given total cm. energy W. In terms of the T J the total helicity amplitude 
for a — > b is given by 

Tx ai x a .^x bl x b2 (WA,0) = 87rW£(2J+ lJe'^X^W^^^^ (3) 

j 

with A a = A ai — A a2 , A;, = A^ — A& 2 , and 9, <fi are the spherical coordinates of the outgoing 
particle b\. 

We apply relation @ to the three coupled channels 7r°d, pn, and 70? at the 7r° production 
threshold, i.e., in the limit W — > W = m n o +rrid, which, of course, implies p n o d — > 0. In this 
limit, therefore, no contribution from c = (ir°d) to the sum occurs due to the vanishing phase 
space factor. Moreover, at threshold we can restrict ourselves to partial waves of angular 
momentum and parity J n = 1~. Thus 7i°d — > pn is described by a single matrix element A, 
with 

Tt d -> XpK = ^(101A P - A n |lA p - \ n )A. (4) 

For 7c? — * TT°d at threshold, two independent matrix elements remain, E n od and M^o^. They 
correspond to the electric dipole (El) and magnetic quadrupole (M2) radiation allowed for 
the 1 + — > 1 _ transition of the hadronic system. One finds 

T l^' d = V7^T, v / 2TTT(l - A d LA 7 |lA 7 - A,) (5 L>1 E^o d + X^M^). (5) 

3 



Finally, for 7c? — > 3 Pi(pn) one obtains 
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£ V2L + 1 (1 - A d LA 7 |lA 7 - A d ) (6 LA E pn + A 7 5 L , 2 M pn ), (6) 

£=1,2 

where E pn and M pn are the corresponding El and M2 matrix elements of the disintegration 
process. The various total cross sections in terms of these matrix elements are 

<r{^d^ P n) = ^^\A\\ (7) 
a(jd - n°d) = (Wo/ + |M^| 2 ) , (8) 

a ( 7 d - 3 Pi(pn); E1+M2) = (|£ pn | 2 + |M p „| 2 ) . (9) 

Taking now a = (70?) and b = (7r°d), relation (0) leads to 

QfmEMWo) = ■j=p pn E pn (Wo)A*(W ), (10) 



and an analogous relation for M n o d and M pn , respectively. Since the lhs of ( ID ) is real, 
relation (|H]) implies that the phases of the complex matrix elements A and E pn are equal. 
Indeed, evaluating (||) with a = (vr c/), 6 = (pn) and a = (7c?), b = (pri) provides us with 

A(W ) = \A(W )\ exp(i83 Pl (W )+imr), (11) 

and 

E pn (W) = \E pn (W)\ exp (t5 3Pl (W) + ikir) , W < W , (12) 

respectively, where Ssp 1 is the nucleon-nucleon scattering phase shift in the 3 Pi channel. Eq. 
(|12|) is simply Watson's theorem applied to deuteron photodisintegration and is valid for 
all energies below the pion production threshold, whereas Eq. ( |TTD is valid for W = W$ 
only. Finally, we mention that taking a = b = (nd) in (0), leads to the constraint for the 
imaginary part of the Tid scattering length in ([!]). 

Using (|7]) and (|), and the detailed balance relation 
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3plo d a(7v d ^ pn) = 4p 2 pn a(pn -> 7r°d), (13) 
our main result fllPP can be rewritten as 



IQfm-E^dl = —^-p pn J a (jd 3 Pi(pn); El) (j(pn -> 7r°d). (14) 
V27T y p n od 

At this level, there is, however, no way to fix the sign. In order to get a numerical value, we 
take for the hadronic cross section the experimental result given by Hutcheon et al. J7|, 

lim 2 — ^ a{mi -> vr°rf) = 184 ± 5 ± 13 //&. (15) 

The partial cross section 7<i — > 3 Pi(pn) is at present not available although it could in prin- 
ciple be obtained from a multipole analysis. There are, however, reliable theoretical models 
available which reproduce all deuteron photodisintegration data in this energy region |§. 
At this energy, the El matrix element is mainly given by 7r-exchange current contributions, 
which can be calculated in a largely model-independent way by taking advantage of gauge- 
invariance constraints (Siegert's theorem). The underlying nucleon-nucleon interaction and 
also all model-dependent transverse electromagnetic currents, like the A(1232) excitation 
current, have little effect on this matrix element. We take the value 



a 



(jd -> 3 Pi(pn); El) = 10.5//6 (16) 



from an updated version of the model of Ref. 0. In order to relate our result to the 
calculations of Beane et al. Q3] , we switch to their normalization (and notation) of the electric 
dipole amplitude, Ed = ^0^/4, and obtain 

\SmE d \ = 0.22 x 10- 3 /m w +. (17) 

The cross section a (jd — > 3 Pi(pn); M2) is more than two orders of magnitudes smaller than 
(|TJ) and leads to |3m M d \ = 0.018 x 10- 3 /m n +, where M d = M^d/A. Nevertheless, the 
role of the M2 transition in the coherent production at threshold deserves a more detailed 
investigation which will be presented elsewhere. The main reasons are: (i) it provides 
an additional possibility to test theoretical predictions (an experimental separation of El 



and M2 requires a polarized deuteron target), and (ii) the relative importance of the M2 
transition grows with the momentum transfer in the electroproduction process. 
The result ([T7|) has to be compared with the value calculated in [|j], 

Ef T = 0.38Ef + n - 2.6 x lO^/m^ = -1.8 x 10- 3 /m w+ , (18) 



where the latter value is obtained taking the %PT prediction of [ TOJ for the neutron electric 
dipole amplitude, E^ 1 = 2.13 x 10~ 3 /m w +. Thus the threshold cross section itself is affected 
by less than 2% only by this imaginary part (|17D . However, there is no reason at all to assume 
that the contribution of the absorptive process to the real part of the amplitude is much 
smaller in magnitude than the imaginary part. Unfortunately, unitarity does not allow to 
estimate it. It can only be calculated within a model which to our knowledge has not yet 
been done. For the moment, in order to get a rough idea, one may look into the in many 
respects analogous situation for ird elastic scattering. There, three-body calculations suggest 
for the absorptive contribution to the scattering length, 9fte ~ —Qm a w d (see || for an 
overview). Assuming, therefore ?R.e E°[ bs = 0.22 x 10 _3 /m 7r +, one has to conclude that the 
neglect of the absorption process in an analysis of the 7c? — > 7r°d cross section based on ( |T8"D 
would lead to a systematic error of the order of SEq'^J 1 = ±0.6 x 10~ 3 /m 7r + for the neutron 
amplitude. This rough argument at least demonstrates that a calculation of the absorptive 
contribution is necessary before definite conclusions on the neutron amplitude can be drawn. 
One way to do this is to combine conventional models for deuteron photodisintegration || 



with those for the pionic disintegration of the deuteron JTT]. 

As second point we would like to address the question what is the role of this absorptive 
process in the framework of heavy baryon ;\TT. The xPT treatment of the process 7c/ — > 7i°d 
is sketched in Fig. [1]. The non-absorptive contribution, Fig. |l](a), is based on the two-nucleon 
irreducible kernel K^ n for the jpn — > 7T°pn subprocess. It is obtained from the effective chiral 
Lagrangian, and then sandwiched between deuteron wave functions ipd- K J7T sums all time- 
ordered diagrams which do not contain pure two-nucleon intermediate states. It has been 
calculated in Q up to and including all terms of order v < of the expansion in terms of 



powers of small momenta (Q/A) u where typically Q ~ m n , A ~ m^. 

The proper treatment of the absorptive contribution is shown in Fig. |l](b). Formally it 
is given by 

T xPT/a bs = KvGqKi ^ + Kn G T NN (W + ie)G K 1 (19) 

Here, the input from xPT are the two-nucleon irreducible kernels and K n for '-/pn — > 
and — > ir°pn, respectively. These are linked by either the free two-nucleon propagator, 

Go = (Wo - 2m N - f/m N + ie) , (20) 

or the propagation via G T NN G which includes the nucleon-nucleon interaction through 
the full off-shell T-matrix in the 3 Pi partial wave. T^n (and ipd) has to be calculated by 
solving the Lippmann-Schwinger equation for a two-nucleon potential. Ideally, the potential 
is thereby also taken from xPT as the sum of all irreducible NN — > NN diagrams. 

Now the question arises, where the absorptive contribution fits into the %PT power 
counting scheme. We will answer it by means of Fig. The time-ordered graph of Fig. 0(a) 
shows an absorptive contribution which is contained in the first term of flrS|). The diagram 
of Fig. 0(b) is a part of K ln build up from the same interaction vertices as (a). According 



to Weinberg's power counting rules |12Hl^| , irreducible graphs with N nucleons, C separate 



connected pieces, L loops, and Vi vertices of type % contribute to order Q v with v given by 

v = A-N -2C + 2L + ^y i A i , (21) 

% 

where the index A« is bounded by chiral invariance. For pure hadronic vertices one has 
Aj > 0, while for vertices with one photon Aj > — 1. For the graph of Fig. Q(b) with N = 2, 
C = L = 1 this leads to 1/2(6) > 1- Consequently these type of contributions could be ignored 
by Beane et al. in order to be consistent up to and including v < 0. 

However, as has been stressed by Weinberg, graphs with intermediate states containing 
only nucleons violate the simple power counting rules, because of the small (nearly infrared- 
divergent) energy denominators associated with the propagation of these states, see Eq. 



(pO|). Fig. 0(a) belongs to this class of reducible graphs. In time-ordered perturbation 
theory, energy denominators of states with at least one pion are of order Q, while those with 
nucleons only are of order Q 2 /rriN, when one assumes that all momenta are of order Q. This 
leads to count the order of Fig. 0(a) as u 2 ( a ) = v<i{b) — 1 > 0. 

It could be objected, that the momenta of the intermediate particles in Fig. 0(a) are 
not of order Q but rather of order P ~ ^Jm^m^ which implies that the two-nucleon energy 
denominator is of order Q~ x rather than Q~ 2 . Actually this is true for the imaginary part of 
the diagram which just arises from the kinematical situation where the nucleon momenta are 
equal to the on-shell momentum p pn which is given by V 2 P nl m N = m ^ applying nonrelativistic 
kinematics. In such a situation also the two intermediate pion momenta must be of order 
P in order to change the relative nucleon momenta of order Q, provided by the deuteron 
wave function. However, as will be seen soon, counting all powers of P in the graph (a) of 
Fig. Q will not change the above conclusion that it contributes already to the order Q°. The 
necessity to consider a modified power counting (due to the kinematics of the reaction) was 
first noted by Cohen et al. flT5|, studying the reaction pp — > 7i°pp near threshold. In order 



to get the worst case, we assume from now on that all vertices in Fig. 0(a) arise always 
from the leading terms of the chiral Lagrangian. Following the steps in |T2|] one counts for 
Fig. 0(a): P 3 from three derivative couplings, P~ 2 from the two energy denominators of 
the states containing pions, P~ 2 from four factors 1 / y/2E n , a factor P 3 from the integral 
over the loop three-momentum, and finally P~ 2 from the two-nucleon energy denominator. 
Altogether, the graph (a) of Fig. gives a contribution of the order P° or equivalent Q°. 

Thus we have to conclude that a complete calculation including all terms of order Q u 
with v < requires the inclusion of the absorptive process. Indeed, already the imaginary 
contribution of the absorption process, ^smE^ in ([H]), turns out to be of the same order 
of magnitude as the three-body contribution of order v = calculated by Beane et al., 
£* M = -0.25 x l(r 3 /m^+ in the notation of §. 

In summary, the contribution of the two-step process 'yd — > pn — > ir°d to the imaginary 
part of the electric dipole (and magnetic quadrupole) amplitude for coherent pion photo- 

8 



production on the deuteron has been calculated utilizing unitarity constraints. The result 
shows that this absorptive process cannot be neglected in the extraction of the elementary 
neutron amplitude. It has been shown that a consistent calculation for 7c? — > -n°d 

beyond next-to-leading order requires indeed the inclusion of the absorptive process. 

I thank H. Arenhovel for useful discussions. This work was supported by the Deutsche 
Forschungsgemeinschaft (SFB 201). 
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(a) (b) 
FIG. 1. Schematical representation of the 7a! — > 7r°a! amplitude at threshold in the xPT 

framework: the non-absorptive part (a) is based on the irreducible jpn — ► 7r°p"- kernel -ftT 77r , and 

the absorptive part (b) combines the irreducible kernels K 7 and K w for 7j»n — > pn and — ► n pn, 

respectively, ipd is the deuteron wave function and T/vat the nucleon-nucleon scattering matrix in 

the 3 P\ partial wave. 





(a) (b) 
FIG. 2. Two diagrams to illustrate our power counting arguments. The time-ordered diagram 

(a) is part of the absorptive process, whereas diagram (b) although built up from the same vertices 

contributes to the non-absorptive part due to a different time-ordering. Solid, dashed, and wavy 

lines represent nucleons, pions, and photons, respectively. 
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